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ABSTRACT 

We present an occultation of the newly discovered hot Jupiter system WASP-19, ob- 
served with the HAWK-I instrument on the VLT, in order to measure thermal emission 
from the planet's dayside at ^2\xm. The light curve was analysed using a Markov- 
Chain Monte-Carlo method to find the eclipse depth and the central transit time. 
The transit depth was found to be 0.366±0.072%, corresponding to a brightness tem- 
perature of 2540±180K. This is significantly higher than the calculated (zero-albedo) 
equilibrium temperature, and indicates that the planet shows poor redistribution of 
heat to the night side, consistent with models of highly irradiated planets. Further ob- 
servations are needed to confirm the existence of a temperature inversion, and possibly 
molecular emission lines. The central eclipse time was found to be consistent with a 
circular orbit. 

Key words: methods: data analysis, stars: individual (WASP-19), planetary systems, 
techniques: photometric 



1 INTRODUCTION 

The vast majority of exoplanets are discovered due to effects 
they have on the observed orbits and brightness of their host 
stars. Transiting planets are detected when the planet passes 
in front of its host star, causing the amount of light we detect 
to decrease periodically. From the shape of the transit light 
curve, coupled with radial velocity measurements, we can 
determine the radius and mass of the planet, and infer its 
composition. 

During occultation, when the planet passes behind its 
host star, we have the opportunity to directly measure the 
flux emitted by the planet. This allows detection of thermal 
emission from the planet's dayside and reflected starlight, 
and in the near-infrared (where the flux from the planet is 
dominated by thermal emission), yields a measurement of 
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the planet's temperature. Occultations of transiting planets 
were first detected from space using the Spitzer space tele- 
scope ( Charbonneau et al.|2005 Deming et al.|2005 l, which 
has continued to measure thermal emission at wavelengths 
3.6 ii.m. Ground based infrared photometry should pro- 
vide complementary information, providing measurements 
of thermal emission at shorter wavelengths. This is the re- 
gion where the flux from a typical hot Jupiter is expected 
to reach its peak, and also where molecular absorption (or 
emission) bands in the atmospheres play an important role 
(see e.g. 



Barman et al. 
et al.|2008| |. 



2005 Marley et al. 2007 Fortney 



However, ground based infrared photometry of occulta- 
tions has proven extremely difficult to date, due to the very 
high signal-to-noise requirements and the relatively poor 
stability of IR detectors in comparison to optical CCDs. 



Following some tentative detections (e.g. Snellen & Covino 



20071, the first ground based detections of thermal emission 



were recently published by de Mooij fc Snellen ( 2009 , TrES- 



3b in K-band),^ Sing fc L6pez-Morales| ( |20(j9| OGLE-TR-56b 

2]iin). 



in z' band) and Gillon et al 
Here 



(2009 



CoRoT-lb at 
we present observations of an occultation of 
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the newly discovered transiting exoplanet system WASP- 19 
( Hebbet al.|2010[ ). WASP- 19 is a G-dwarf star with an eflfec- 
tive temperature of 5500±100K, with a 1.15±0.08 Jupiter- 
mass planet in a very short period orbit of 0.7888399 ± 
0.0000008 days. The planet is therefore expected to be ex- 
tremely hot due to the proximity to its host, and coupled 
with its large radius of 1.31 ± 0.06 Jupiter-radii is expected 
to produce an occultation in the K-band ~ 0.2 — 0.4%. It 
was observed using the HAWK-I (High Acuity Wide-field 
K-band Imager) instrument ( Pirard et al.|2004 Casali et al 



2006 Kissler-Patig et al. 20081, a wide-field, near-infrared 



camera mounted on UT4 of the VLT, which was recently 
used by Gillon et al. (2009) to detect the occultation of 



CoRoT-lb. Simultaneously with the original submission of 



the present manuscript, Anderson et al. (20101 announced 



an independent detection of the secondary eclipse of WASP- 
19b in the H-band, also with HAWK-I, finding an eclipse 
depth of 0.259t°Q°Qll%. In Section [2j we describe the obser- 
vations and data reduction. Sections [S] and [4] present our 
results and analysis, and finally in Section [5] we summarise 
and discuss our findings. 



2 HAWK-I OBSERVATIONS AND DATA 
REDUCTION 

An occultation of the exoplanet system WASP-19 was ob- 
served using HAWK-I. The observations were carried out in 
service mode on the night of March 30 2009, and lasted for 
~4.6 hours, starting ~1.5 hours before the expected ingress, 
and ending ~1.6 hours after expected egress (assuming the 
planet is in a circular orbit). HAWK-I consists of four Hawaii 
2RG 2048x2048 pixel detectors, arranged in a mosaic giving 
a total field of view of 7.5 arcminutes squared, with a pixel 
scale of 0.106 arcseconds. However, in order to obtain very 
precise relative photometry needed to detect a secondary 
transit, we found (similar to Gillon et al.||2009l ) that only 
reference stars that fell on the same chip (quadrant 2) as 
WASP-19 could be used, reducing the effective field of view. 
An image of the HAWK-I field (quadrant 2 only) is shown 
in Fig.[fl 

Observations were taken using the NB2090 filter, a 
narrow filter centred at 2.095 (xm, in the region where the 
thermal emission of hot Jupiters (~ 1000 — 1800 K) is ex- 
pected to peak. Using a narrow filter should also help re- 
duce differential extinction. This still provides high enough 
signal-to-noise to detect the occultation assuming observa- 
tions are photon limited, in which case we would obtain sub- 
millimagnitude photometry per 10 second exposure]^ How- 
ever, from experience we know observations of this type are 
limited by systematic noise. The telescope was defocussed 
slightly to keep the stars well within the linearity range of 
the detectors, and to spread the stellar profile over more 
pixels. 

A total of 185 images were acquired, each consisting 
of 6 X 10 second exposures. A dithering pattern was cre- 
ated at random within 30 arsceconds of the first image us- 
ing the AutoJitter template. A basic calibration including 



^ According to the VLT/HAWK-I exposure time calculator 
http://www.eso.org/observing/etc/ with K = 10.5. 



Figure 1. HAWK-I image of quadrant 2 showing the position 
of WASP-19 (box) and the five comparison stars used for relative 
photometry (circles). The box around WASP-19 shows the size of 
the dithering box used by the AutoJitter template. 




Figure 2. Close-ups of WASP-19 in the HAWK-I images showing 
how the PSFs elongated during the course of the observations. 
The left image shows a typical PSF during the first half of the 
night with a round shape. The right image shows a typical PSF 
towards the end of the night, with an elongated shape. 



bias subtraction and fiat-fielding was performed on each im- 
age using the EsoRex HAWK-I pipeline. Sky-subtraction 
was also performed using EsoRex. For each image, a sky 
background frame was constructed from a running median 
of the object frame and the nine previous and nine subse- 
quent frames. This process involved two passes, the second 
of which masked objects from the background calculation. 
The stars had a typical FWHM of 0.5-0.6 arcseconds, but 
around half way through the observations the PSFs of the 
stars deteriorated and began to elongate, as shown in Fig.[2j 
probably related to the defocussing of the telescope. The oc- 
cultation signal could still be recovered from the photometry, 
providing large enough apertures were used. However, this 
is likely responsible for the large gradient seen in the raw 
light curve. 

Aperture photometry was performed on WASP-19 and 
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Figure 3. Raw VLT/HAWK-I light curve of tiie secondary tran- 
sit of WASP-19. Tlie dashed-dotted lines show the expected start 
and end of transit, assuming the planet is in a circular orbit. A 
transit is clearly visible at the expected time and duration. The 
solid and dashed lines show the best-fit models using the linear 
and quadratic normalisation functions, respectively. The first and 
last 20 images were not used in the fitting process, and the cor- 
responding data points are shown in grey. The error bars shown 
were re-scaled so that the recuded of the best-fit model (using 
the linear normalisation) is equal to one. 



five comparison stars using PyraQand tlie DAOPHOT pack- 
age. Tlie positions of WASP-19 and tiie five comparison stars 
used are marked in Fig. [l] Several brighter comparison stars 
were ruled out either because they were over the linearity 
limit of the chip, or because they feU too near the edge of the 
field on some images due to the dithering. An aperture of 12 
pixel radius was used for each star, and the remaining sky 
was estimated using an annulus of width 5 pixels beginning 
at 20 pixels from the stars' centres. The flux of WASP-19 
was then divided by the sum of the flux from the flve com- 
parison stars to obtain the raw light curve, shown in Fig. [3] 
Initial estimates of the photometric errors were calculated 
using the aperture electron flux, sky and read noise. 

Clearly, a trend can be seen in the light curve, likely 
caused by the PSFs elongating during the observations. This 
was removed by normalising the light curve by fitting a lin- 
ear or quadratic function of time to the out-of-transit data, 
to set the out-of-transit flux equal to 1. The best-fit mod- 
els using the linear and quadratic functions are shown in 
Fig. [3] although the linear function is preferred and used to 
determine the light curve parameters (see Sect.jsjl. The nor- 
malisation of the light curves for both functions was strongly 
affected by outlying points at the start and end of the obser- 
vations. As the sky background was estimated from 9 frames 
at either side of each image, the sky subtraction was not es- 
timated accurately. We decided to remove the first and last 
20 data points (~ 30 minutes) from the light curve anal- 
ysis (shown in grey in Fig. [3|, which still left substantial 
out-of-transit data to constrain the normalisation function. 



Pyraf is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 
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Figure 4. Top:VLT/HAWK-I light curve of the secondary tran- 
sit of WASP-19 normalised with a linear function of time, with the 
best-fit model over-plotted and with the first and last 20 images 
removed. Bottom: Residuals from the best-fit model. 



3 LIGHT CURVE MODELLING AND 
ANALYSIS 

The HAWK-I light curve was fitted using an MCMC 
(Markov-Chain Monte-Carlo) routine, which is a method 
used to explore the multi-dimensional parameter space of a 
model fit efficiently, allowing a determination of the joint 
posterior probability distribution for the parameters (see 
e.g., ITegrnark et al.| pOOil |Holman et al.| [20061 |Collier| 



Cameron et al. 



2007 



tion of MCMC uses the x 
curve given by 



Winn et al 
•2 



20081. Our implementa- 



fitting statistic for a model light 



JV 

i=i 



(/obs,j /calc.j) 



where /obsj is the flux observed at time j, aj is the corre- 
sponding uncertainty and /caic.j is the flux calculated from 
the model for time j. 

The model flux /caic was constructed using Kepler's 
laws to determine the normalised separation of the planet 
and star centres as a function of time assuming a circular 



orbit. The analytic equations of Mandel & Agol ( 2002 1 were 
then used to calculate the flux from the normalised separa- 
tion assuming no limb darkening on the planet's surface (as 
this will only affect the ingress and egress which the HAWK- 
I data is unable to constrain). The system parameters were 
held fixed at the values determined by |Hebb et al.| ( [2010[ ) to 
determine the shape of the transit. The depth D and width 
W of the transit could then be varied by simply scaling the 
transit shape in both fiux and time axes. 

A long chain of model parameter sets is created by 
adding small Gaussian perturbations to the previous ac- 
cepted set. At each step we apply the Metropolis-Hastings 
rule to decide whether to accept the new parameter set. A 
model that produces a lower than the previous model 
is always accepted. A model that produces a higher is 
accepted with probability exp(— Ax^/2). This maps out a 
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probability distribution which can be used to estimate best- 
fit values and uncertainties for each parameter. More details 
on the MCMC procedure used in this work can be found in 



Table 1. Parameters and la uncertainties as derived from the 
MCMC fits, and the calculated brightness temperature. 



Gibson et al] ( ^2008. 2010j ). 

The light curve was fitted for D and the central transit 
time Tq. We did not fit the transit width W by allowing it 
to vary as a free parameter, but to account for uncertainties 
that may propagate to D and Tq, we allowed the transit 
width to vary within a Gaussian prior by adding a term to 
the statistic so that 



2 _ (.fobs.j — /calc.j)^ 



where aw was set to the fractional error in the transit du- 
ration from iHebb et al. (20101. This led to no significant 



changes to the parameters and uncertainties determined, 
but nonetheless we used it in our final analysis for com- 
pleteness. We allowed the stellar flux to vary either linearly 
or quadratically as a function of time to normalise the light 
curves, using a further 2 or 3 normalisation parameters. In 
order to account for possible correlations between these nor- 
malisation parameters and the transit parameters, the nor- 
malisation parameters were allowed to vary freely during the 
fitting process. 

An initial MCMC analysis of length 20 000 was used to 
estimate the jump functions for D,To,W and the normalisa- 
tion parameters. The photometric errors aj were re-scaled 
so that the best-fitting model had a reduced equal to 
one, as there are often further sources of noise not taken 
into account from photon statistics alone. This required the 
errors to be re-scaled by factors of 6.17 and 6.14 when using 
linear and quadratic normalisation functions, respectively, 
showing the quadratic function provides a marginally better 
fit. Such large factors are likely due to variations in pixel-to- 
pixel sensitivity, which are exaggerated by the dithering pat- 
tern required to obtain accurate sky estimates for infrared 
photometry. 

We also examined the residuals for evidence of time- 



correlated noise using the method of Winn et al. (20081, 
where a factor 1) is calculated by analysing the residuals 
from the best-fit model. The noise should drop by ~ 1/ y/~N 
if it is uncorrelated in time, when the data is averaged into 
bins of width A''. However, the noise is often larger by a fac- 
tor /3 indicating an extra time-correlated source of noise. See 



Gibson et al. ( 2008 1 for a more detailed description of this 
procedure. We found relatively low levels of time-correlated 
noise, with the maximum value of /3 ~ 1.08 for both the 
linear and quadratic fits, when the residuals were binned 
into ~ 12 minute intervals. This indicates that the noise is 
probably dominated by systematic errors in the randomised 
dithering pattern that aren't correlated in time. The photo- 
metric errors CTj were then rescaled by P prior to the final 
MCMC analysis. 

The full MCMC analysis was then performed consisting 
of five separate chains each of length 100 000. The first 20% 
of each chain were eliminated to keep the initial conditions 
from influencing the results. The adopted value of each pa- 
rameter was taken to be the modal value of the probability 
distribution, and the uncertainties taken as the limits which 
encompass the central 68.2% of the probability distribution. 
To check that the chains had all converged to the same prob- 



Parameter 


Symbol 


Value 


Unit 


Transit depth 


D 


0.366±0.072 


/o 


Transit centre 


To 


2454921. 66790±0.00190 


HJD 


Brightness temp 


Tb 


2540±180 


K 



abiUty distribution, the Gelman & Rubin Statistic ( [Gelman] 
fc Rubin|[T992| ) was calculated for each of the free parame- 
ters and found to be less than 0.5% from unity for each, a 
good sign of mixing and convergence. 

As using the quadratic normalisation function only pro- 
duces a marginally better fit to the data, it is useful to 
test whether we are 'over-fitting' the data. To check which 
model to use, we evaluated the Bayesian Information Crite- 
rion (BIG, |Schwarz||1978[ ) for the best-fit model using both 
the linear and quadratic normalisation functions. The BIG 
is given by 

BIC^2 =x^ + fcln(n), 

where is calculated for the best-fit model, k is the num- 
ber of parameters in the model fit, and n is the number of 
data points. Increasing the number of free parameters always 
gives a better fit, but the BIG penalises for higher numbers 
of model parameters to test the model is not over-fitting 
the data, and the model that gives the lowest BIG value is 
the preferred model. Calculating this for the both the linear 
and quadratic normalisation, indicates that the quadratic 
normalisation does not improve the fit sufficiently to justify 
its use, and the results from the model using the linear nor- 
malisation function are the adopted measurements. This is 
true when is calculated using the photometric error bars 
both before and after they are re-scaled. A plot of the nor- 
malised light curve and resulting best-fit model using the 
linear normalisation function is shown in Fig. |4] 



4 RESULTS 

Results from the light curve fits are shown in Table [l] The 
transit depth is found to be 0.366±0.072%. As a secondary 
check on the transit depth, the depth was found by combin- 
ing the in-transit and out-of-transit data points of the nor- 
malised light curve and calculating the depth from the differ- 
ence. This gave D = 0.350±0.065%, in agreement with the 
depth derived from the MGMG analysis, but with smaller 
uncertainty as this doesn't take into account uncertainties 
in the central transit time, transit width, and normalisation 
parameters. 

The central transit time was determined to be 
2454921. 66790±0.001 90 HJD. Using the ephemeris of 
Hebb et al. (20101, this corresponds to a phase of 
0.50114±0. 00241, showing that the central eclipse time is 
consistent with a circular orbit. 
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5 SUMMARY AND DISCUSSION 

We have detected an occultation of the extrasolar planet 
system WASP-19 in the NB2090 filter centred at 2.095 \im 
using the HAWK-I instrument. The eclipse depth was mea- 
sured to be 0.366±0.072%, with the transit centre occurring 
at phase 0.50114±0. 00241, consistent with a circular orbit to 
well within 1 cr. [Anderson et al.| ( |2010| ), however, report evi- 
dence for a non-zero eccentricity at the 2.6(t level. The meth- 
ods used by Anderson et al.| ( 2010 1 to measure the phase of 
the eclipse, and to correct for systematics, particularly the 
trends in the light curve which can considerably effect tran- 
sit times (e.g. Gibson et al.|2069 l, differs from that adopted 
in our analysis, and this difference is the most likely source 
of the apparent discrepancy. 

The brightness temperature may be calculated from 
the depth, taking the values and uncertainties for the stel- 
lar temperature and planet-to-star radius ratio from |Hebb| 
et al. (20101. This results in a brightness temperature 



of 2540±180K, considerably larger than the (zero-albedo) 
equilibrium temperature given in Hebb et al. (20101 of 



2009±26K. Anderson et al. (20101 report an eclipse depth 
of 0.2591q q44% in the H-band, corresponding to a bright- 
ness temperature of 2560±130K. This is again considerably 
higher than the equilibrium temperature, and also the equi- 
librium temperature given no redistribution of heat to the 
nightside, which they calculate as ~ 2400 K. Our 2.095 \im 
measurement is consistent with this temperature. 

This indicates that there is poor redistribution of heat 
to the night side of the planet, consistent with the pM 



class of planets in the classification of Fortney et al. ( 2008 1 



They also conclude that planets hot enough to have sig- 
nificant TiO and VO absorption in their upper atmospheres 
would show temperature inversions, and may appear anoma- 
lously bright in the infrared due to molecular emission 
bands, which may help explain the temperature excesses 
here. Other ground based K-band secondary eclipse mea- 
surements of very hot Jupiters reached similar conclusions 
(|de Mooij fc Snellen|[2009} [Gillon et aLl[2009l [Aionso et al.| 
2010). Further observations at infrared wavelengths are re- 



quired to confirm temperature inversions and possibly mea- 
sure molecular emission for WASP-19, which should prove 
an interesting target for future studies. 
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